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Abstract 

Neutral Higgs bosons of the Standard Model (SM) and the Minimal Super- 
symmetric Standard Model (MSSM) were searched for in the data collected in 
1999 by the DELPHI experiment at centre-of-mass energies between 191.6 and 
201.7 GeV with a total integrated luminosity of 228 pb _1 . These analyses, in 
combination with our results at lower energies, set 95% confidence level lower 
mass bounds on the Standard Model Higgs boson (107.3 GeV/c 2 ) and on the 
lightest neutral scalar (85.9 GeV/c 2 ) and neutral pseudoscalar (86.5 GeV/c 2 ) 
Higgs bosons in representative scans of the MSSM parameter space. An ex- 
tended scan of the MSSM parameter space was also performed to test the 
robustness of these limits. 
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1 Introduction 

In the framework of the Standard Model (SM) there is one physical Higgs boson, H, 
which is a neutral CP-even scalar. At LEP2 the main production process is through the 
s-channel, e + e~— > Z* — >HZ, but there are additional t-channel diagrams in the B.uu and 
He + e~ final states, which proceed through W + W~ and ZZ fusion, respectively. With 
the data taken previously up to y/s = 188.7 GeV, DELPHI excluded a SM Higgs boson 
with mass less than 94.6 GeV/c 2 at the 95% confidence level (CL). The other LEP 
collaborations reached similar results J2|. The present analysis concentrates on masses 
between 85 and 115 GeV/c 2 . The results obtained in the same mass range with the 
data taken by DELPHI in the last year of LEP operation and analysed with preliminary 
calibration constants can be found in ||. Although the emphasis is on high masses, the 
analysis described in this paper is also applied to lower masses, down to the bb threshold, 
in order to derive a constraint on the production cross-section of a SM-like Higgs boson 
as a function of its mass. 

In the Minimal Supersymmetric Standard Model (MSSM), the production of the 
lightest scalar Higgs boson, h, proceeds through the same processes as in the SM. The 
results of the search for the SM Higgs boson thus also apply to the h boson. However, 
in the MSSM the production cross-section is reduced with respect to the SM one and 
even vanishes in part of the MSSM parameter space. This model predicts also a CP-odd 
pseudo-scalar, A, that would be produced mostly in the e + e~— > Z* — > hA process at 
LEP2. The MSSM parameters are such that when the single h production is suppressed, 
the associated hA production is enhanced. This channel is thus also considered in this 
paper. Previous 95% CL limits from DELPHI on the masses of h and A of the MSSM 
were 82.6 GeV/c 2 and 84.1 GeV/c 2 respectively J|. The results of the other LEP collab- 
orations are described in Ref. f2j. The present analysis in the hA channel concentrates 
on masses between 80 and 95 GeV/c 2 . 

In the HZ channel, all known decays of the Z boson have been taken into account 
(hadrons, charged leptons and neutrinos) while the analyses have been optimised either 
for decays of the Higgs particle into bb, making use of the expected high branching 
fraction of this mode, or for Higgs boson decays into a pair of r's, which is the second 
main decay channel in the SM and in most of the MSSM parameter space. A dedicated 
search for the Higgs boson invisible decay modes will be reported separately. The hA 
production has been searched for in the two main decay channels, namely the 4b and 
bbr + T~ final states. 

2 Data samples and detector overview 

In 1999 LEP ran at centre-of-mass energies ranging from 191.6 GeV to 201.7 GeV. 
DELPHI recorded 25.9 pb" 1 at 191.6 GeV, 76.9 pb" 1 at 195.6 GeV, 84.3 pb~ x at 
199.6 GeV and 41.1 pb _1 at 201.7 GeV. The requirement of full detector performance 
reduces the luminosities in the He + e~ and Hi/z/ searches by at most 3%. The detector was 
unchanged from the previous data taking period. Ref. [Q provides a short description 
while more details can be found in Ref. [||,^] for the original setup and in Ref. |7] for the 
LEP2 upgrade of the silicon tracking detector. 

Large numbers of background and signal events have been produced by Monte Carlo 
simulation and then passed through the DELPHI detector simulation program || . These 
samples typically correspond to about 100 times the luminosity of the collected data. 
Backgrounds were generated with PYTHIA || for hadronic two-fermion final states (here- 
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after denoted as qq(7)) and with KORALZ |J for leptonic two-fermion final states. The four- 
fermion background, which is a coherent sum of many processes whose main components 
are referred as Z7*, W + W~ and ZZ in the following, was generated with EXCALIBUR |10 



in most of the phase space, but GRC4F JO] and KORALW [12j were used to complete the 
EXCALIBUR samples in the case of very forward electrons or low mass hadronic resonances, 
respectively. TWOGAM and BDK [14" were used for two-photon processes (hereafter de- 
noted as 77) and BHWIDE [[TJJ for Bhabha events in the main acceptance region. 

Signal events were produced using the HZHA [IB] generator. As the enhancement in the 
production cross-section due to W + W~ fusion and to its interference with the HZ process 
is significant in the Huu channel [17|] (for mn around the HZ kinematic threshold), signal 
events in this channel were generated using a version of HZHA modified to include also 
fusion and interference between the HZ and W + W~ fusion diagrams. For the HZ process, 
the H mass was varied from 12 to 115 GeV/c 2 , while for hA the range for the A mass 
was 12 to 95 GeV/c 2 . A step of 5 GeV/c 2 was used above 80 GeV/c 2 , since the analyses 
were optimized at high mass. Wider steps were used at lower masses. Moreover, the hA 
signal events were simulated for three values of tan f3 (the ratio of the vacuum expectation 
values of the two Higgs field doublets of the MSSM) equal to 2, 20 or 50. This fixes the h 
mass, almost equal to wa for tan/3 = 20, 50 and lower than by around 20 GeV/c 2 if 
tan (3 = 2. The h and A widths are lower than 1 GeV/c 2 for tan j3 below 20 and increase 
rapidly to reach several GeV/c 2 at tan/5 = 50, thus going above the experimental mass 
resolution which is typically of 3 GeV/c 2 in the hA channels. 

The HZ simulated samples were classified according to the Higgs and Z boson decay 
modes. For He + e~, H/x+yU - and Hvv the natural SM mix of H decay modes into fermions 
was permitted. As final states with hadrons and two r's benefit from a dedicated analysis, 
the rr decay mode was removed in the Hqq channel simulations and we generated 
separately the two HZ channels involving r leptons for which one of the bosons is forced 
to decay to r's and the other hadronically. Finally, the hA simulations cover final states 
involving either four b quarks or two b quarks and two r's, irrespective of which Higgs 
boson decays into r's. Efficiencies are defined relative to these states. The sizes of these 
samples vary from 2000 to 3000 events and they were produced at the four centre-of-mass 
energies. 

Although the above signal simulations cover most of the expected final states in the 
SM and MSSM, they were completed by two additional sets at 199.6 GeV. We generated 
hA samples with large mass differences between the h and A bosons, as expected when 
scanning the MSSM parameter space more widely than in the representative scans, and 
hZ samples with h — > AA, as expected in restricted regions of the MSSM parameter 
space. In these two sets, the A (h) mass was varied from 12 GeV/c 2 (50 GeV/c 2 ) up 
to the kinematic limit and only the main decays were simulated; the hA samples were 
restricted to four b final states, and the (h — > AA)Z samples to hadronic decays of the 
Z boson and either four b or four c quarks from the A pair. The results obtained from 
these samples were assumed to be valid also at the three other centre-of-mass energies. 



3 Features common to all analyses 

3.1 Particle selection 

In all analyses, charged particles are selected if their momentum is greater than 
100 MeV/c and if they originate from the interaction region (within 4 cm in the trans- 
verse plane and within 4 cm / sin 6 along the beam direction, where 9 is the particle polar 



3 



angle). Neutral particles are defined either as energy clusters in the calorimeters not as- 
sociated to charged particle tracks, or as reconstructed vertices of photon conversions, 
interactions of neutral hadrons or decays of neutral particles in the tracking volume. All 
neutral clusters of energy greater than 200 or 300 MeV (depending on the calorimeter) 
are used. The 7r =1= mass is used for all charged particles except identified leptons, while 
zero mass is used for electromagnetic clusters and the K° mass is assigned to neutral 
hadronic clusters. 



3.2 b-quark identification 



The method of separation of b quarks from other flavours is described in ||18|| , where 
the various differences between B-hadrons and other particles are accumulated in a single 
variable, hereafter denoted x^ for an event and x l h for jet i. A major input to this 
combined variable is the probability that all tracks with a positive lifetime-signed 
impact parameter^ in the jet lead to a product of track significances as large as that 
observed, if these tracks do originate from the interaction point. A low value of this 
probability is a signature for a B-hadron. The likelihood ratio technique is then used to 
construct x l h by combining P^ with information from any secondary vertex found in the 
jet (the mass computed from the particles assigned to the secondary vertex, the rapidity 
of those particles, and the fraction of the jet momentum carried by them) and with the 
transverse momentum (with respect to the jet axis) of any lepton belonging to the jet. 
The event variable, linear combination of the jet variables. Increasing values of 

Xb (or x\) correspond to increasingly 'b-like' events (or jets). 

The procedure is calibrated on events recorded in the same experimental conditions at 



the Z resonance. The performance of the combined b-tagging is described in Ref. [19] and 



that of the impact parameter tagging alone in Ref. p0| . The overall performance of the 
combined b-tagging for 1999 Z data is illustrated in Fig. [L| Data agree with simulation 
to better than 5% in the whole range of cut values. 



3.3 Constrained fits 



In most channels a constrained fit 21 



is performed to reconstruct the Higgs boson 
mass, and often to reject background processes as well. In order to allow the removal 
of most of the events involving radiative return to the Z, an algorithm has been devel- 
oped p2] in order to estimate the effective energy of the e + e~ collision. This algorithm 
makes use of a three-constraint kinematic fit in order to test the presence of an initial 
state photon along one of the beam directions and hence lost in the beam pipe. This 
effective centre-of-mass energy is called y/s' throughout this paper. 



3.4 Confidence level definitions and calculations 

The confidence level definitions rely on a test-statistic built with the likelihood ratio 
technique p3| . Let Q be the ratio of the likelihood of the observed candidates assuming 



signal plus background to that found using the background-only hypothesis. Q classifies 
the result of an observation between the background-like and signal plus background-like 
situations. We then define the confidence level for the background hypothesis, CLb, as 
the probability, in background-only experiments, to obtain equal or smaller values of 
Q (that is more background-like results) than that observed. Similarly, the confidence 

Throughout the paper, all impact parameters are defined with respect to the reconstructed primary vertex. 
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level for the signal plus background hypothesis, CL s+ b, is the probability, in signal plus 
background experiments, to obtain more background-like results than those observed. 
The pseudo-confidence level for the signal hypothesis, CL S , is conservatively defined as 
the ratio of these two probabilities, CL s+ b/CLb. CL S measures the confidence with which 
the signal hypothesis can be rejected and must fall below 5% for an exclusion confidence 
of 95%. More technical details about how the confidence levels are calculated or how 
uncertainties are taken into account can be found in [Q. 

In the definition of the test-statistic Q, two-dimensional discriminant information is 
used in all channels, as in our previous publication The first variable is the recon- 
structed Higgs boson mass (or the sum of the reconstructed h and A masses in the hA 
channels), the second one is channel-dependent, as specified in the following sections. In 
order to make full use of the information contained in the second variable, the final se- 
lections are loose: the method used for deriving the confidence levels ensures that adding 
regions of lower signal and higher background can only enhance the performance relative 
to a tighter selection, provided the systematic errors are small. 

The distributions are represented as two-dimensional histograms which are derived 
from the simulation samples. These distributions are then smoothed using a two- 
dimensional kernel, which is essentially Gaussian but with a small component of a longer 
tail. The width of the kernel varies from point to point, such that the statistical error on 
the estimated background is never more than 30%. Finally the distribution is reweighted 
so that when projected onto either axis it has the same distribution as would have been 
observed if the smoothing had been only in one dimension. This makes better use of 
the simulation statistics if there are features which are essentially one dimensional, such 
as mass peaks. A check for residual statistical fluctuations was made by dividing the 
simulation into sub-samples, and comparing the expected results; no significant effects 
were observed. 

4 Higgs boson searches in events with jets and elec- 
trons 

The analysis is based upon the same electron identification algorithm and discriminant 
variables as in and is briefly described in the following. The preselection requires at 
least 8 charged particles, a total energy above 0.12-y/J and at least a pair of loose electron 
candidates of energies above 10 GeV and impact parameters below 2 mm (1 cm) in the 
transverse plane (along the beam direction). The Bhabha veto and the modified selections 
allowing for the tau decays of the Higgs boson are as described in JT]] . To reduce the Z7* 
and qq(7) backgrounds, the sum of the di-electron and hadronic system masses must be 
above 50 GeV/c 2 , while the missing momentum is required to be below 50 GeV/c if its 
direction is within 10° to the beam axis. The jet reconstruction and selection proceed as 
in@. 

After this preselection, each pair of electron candidates with opposite charges is sub- 
mitted to further cuts. The electron identification is first tightened, allowing at most 
one electron candidate in the insensitive regions of the calorimeter. The two electrons 
are required to have energies above 20 GeV and 15 GeV. Electron isolation angles with 
respect to the closest jet are required to be more than 20° for the more isolated electron 
and more than 8° for the other one. A five-constraint kinematic fit is performed to test 
the compatibility of the e + e~ invariant mass with the Z mass; the fit imposes energy 
and momentum conservation and takes into account the Breit-Wigner shape of the Z 
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resonance |L]. Events with a fit probability below 1CT 8 are rejected. As the search is 
restricted to high mass Higgs bosons produced in association with a Z particle, the sum 
of the fitted masses of the electron pair and of the hadronic system is required to be above 
150 GeV/c 2 and their difference in the range from -100 GeV/c 2 to 50 GeV/c 2 . The fitted 
hadronic mass and the b-tagging variable Xb are used in the two-dimensional calculation 
of the confidence levels. 

The effect of the selections on data and simulated samples are detailed in Tables [I] 
to while the efficiencies at the end of the analysis are reported in Tables |5| and |5] as a 
function of mn . The agreement between data and simulation at the preselection level is 
illustrated in Fig. § which shows the distributions of the electron energies, the fitted mass 
of the jet system and the isolation angle of the more isolated electron candidate. At the 
end of the analysis, 11 events are selected in the data for a total expected background of 
11.5 ± 0.2(stat.) events coming mainly from the e + e~qq process. 

The systematic uncertainties on background and efficiency estimates are mainly due to 
the imperfect simulation of the detector response and were estimated as described in 
The relative error on the efficiencies is ±3% while that on the background estimates at 
each centre-of-mass energy is ±7%. 

5 Higgs boson searches in events with jets and muons 

The analysis follows that published in , with slight modifications in the preselection 
to adapt to somewhat different beam conditions in 1999. The preselection requires at least 
9 charged particles with two of them in the central part of the detector (40 ° < 9 < 140 °) 
and at least two high quality tracks of particles with a transverse momentum greater 
than 5 GeV/c. For high quality tracks, impact parameters less than 100 /im in the 
transverse plane and less than 500 /im along the beam direction are required. The rest 
of the preselection is unchanged and requires at least two particles of opposite charges 
and momenta greater than 15 GeV/c. 

The rest of the analysis is based upon the same muon identification algorithm and 
discriminant variables as in [|J, but the selection criteria have been re-optimised [Q. As 
a result, the level of muon identification corresponds now to an efficiency of 88% and a 
misidentification probability of 8.8% per pair of muon candidates. At least two muons are 
required with opposite charges, an opening angle larger than 10°, and momenta greater 
than 34 GeV/c and 21 GeV/c. The jet reconstruction and selection proceed as in Q. 
Finally, the angle with respect to the closest jet axis must be greater than 9° for the 
more isolated muon and greater than 7° for the other one. A five-constraint kinematic 
fit taking into account energy and momentum conservation and the Breit-Wigner shape 
of the Z resonance is then performed to test the compatibility of the di-muon mass with 
the Z mass in a window of ± 30 GeV around the Z pole. Events are kept only if the 
fit converges in this mass window. As in the electron channel, the fitted mass of the 
hadronic system and the b-tagging variable Xb are chosen as the discriminant variables 
for the two-dimensional calculation of the confidence levels. 

The effect of the selections on data and simulated samples are detailed in Tables [I| 
to [|, while the efficiencies at the end of the analysis are reported in Tables [5| and |5| as 
a function of m H . The agreement of simulation with data is quite good, as illustrated 
at preselection level in Fig. which shows the multiplicity of the charged particles, the 
momentum of the higher-momentum particle in any preselected pair, the isolation angle 
of the more isolated particle in any preselected pair and the b-tagging variable x^. At the 
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end of the analysis, 8 events are selected in data in agreement with the total expected 
background of 9.4 ± 0.1 (stat.) events coming mainly from the qq process. 

The imperfect simulation of the detector response leads to systematic errors in back- 
ground and efficiency evaluation. As explained in M, each of the momentum and angular 
cuts was varied in a range given by the difference between the mean values of the simu- 
lated and real data distributions of the corresponding variable at preselection level. The 
muon pair identification level, which is a discrete variable, was modified randomly with a 
probability of 5%, corresponding to the maximum difference observed in muon identifica- 
tion results when comparing data with simulation. This is the main source of systematic 
uncertainty in this channel. As a result, a relative error of ±2% is quoted for the efficien- 
cies, independent of mn , while the relative error on the expected backgrounds at each 
centre-of-mass energy is ±3%. 

6 Higgs boson searches in events with jets and taus 

Three channels are covered by these searches, two for the HZ channel, depending on 
which boson decays into t + t~, and one for the hA channel. The analysis, identical to 
that described in |IJ, selects hadronic events by requiring at least ten charged particles, 
a total reconstructed energy greater than 0.4-^/i, a reconstructed charged energy above 
0.2^ and v^s 7 greater than 120 GeV. 

A search for r lepton candidates is then performed using a likelihood ratio technique. 
Single charged particles are preselected if they are isolated from all other charged particles 
by more than 10°, if their momentum is above 2 GeV/c and if all neutral particles 
in a 10° cone around their direction make an invariant mass below 2 GeV/c 2 . The 
likelihood variable is calculated for the preselected particles using distributions of the 
particle momentum, of its isolation angle and of the probability that it comes from the 
primary vertex. As an illustration of the agreement between data and simulation at this 
level of the analysis, Fig. [|a shows the distribution of the isolation angle of the preselected 
charged particle with the highest r likelihood variable in the event. Pairs of r candidates 
are then selected requiring opposite charges, an opening angle greater than 90° and a 
product of the r likelihood variables above 0.45. If more than one pair is selected, only 
the pair with the highest product is kept. The distribution of the highest product of two 
r likelihood variables in the event is given in Fig. f|b. The discrimination between the 
Higgs signal and the SM background is clearly visible. Moreover, the percentage of r 
pairs correctly identified is over 90% in simulated Higgs events. 

Two slim jets are then reconstructed with all neutral particles inside a 10° cone around 
the directions of the r candidates. The rest of the event is forced into two jets using the 
DURHAM algorithm. The slim jets are required to be in the 20°< 9 T < 160° polar 
angle region to reduce the Ze + e~ background, while the hadronic di-jet invariant mass is 
required to be between 20 and 110 GeV/c 2 in order to reduce the qq(7) and Z7* back- 
grounds. The jet energies and masses are then rescaled, imposing energy and momentum 
conservation, to give a better estimate of the masses of both di-jets (r + r _ and qq ), 
that are required to have a rescaled mass above 20 GeV/c 2 , and below y^i to discard 
unphysical solutions of the rescaling procedure. Each hadronic jet must have a rescaling 
factor in the range 0.4 to 1.5. 

The remaining background comes from genuine £ + £~qq events. In order to reject the 
e + e~qq and qq backgrounds the measured mass of the leptonic system is required 
to be between 10 and 80 GeV/c 2 and its electromagnetic energy to be below 60 GeV (see 
Fig. |c). This terminates the selection procedure. The effect of the selections on data 
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and simulated samples are detailed in Tables p] to |4], while the efficiencies at the end of 
the analysis in the three r + r~ qq channels are reported in Tables |5|, [5] and |8] as a function 
of the Higgs boson masses. At the end of the analysis, 6 events are selected in data for 
a total expected background of 6.9 ± 0.2(stat.) events coming mainly from the r + r~qq 
and rz/q'q processes. 

Systematic uncertainties from the imperfect modelling of the detector response were 
estimated by moving each selection cut according to the resolution in the corresponding 
variable. The main contributions arise from the cuts on the t + t~ invariant mass and 
electromagnetic energy. The total relative systematic uncertainties amount to ±6% on 
signal efficiencies and ±11% on the background estimates at each centre-of-mass energy. 

The two-dimensional calculation of the confidence levels uses the reconstructed mass 
given by the sum of the t + t~ and qq di-jet masses after rescaling and a likelihood variable 
built from the distributions of the rescaling factors of the r jets, the r momenta and the 
global b-tagging variable, x\ } . The distribution of this likelihood variable at the end of 
the analysis is shown in Fig. |]d to illustrate the discrimination between the Higgs signal 
and the SM backgrounds. Since the three possible r + r~qq signals are covered by the 
same analysis, the three channels cannot be considered as independent in the confidence 
level computation. For this computation, they are combined into one global t + t~ qq 
channel: at each test point, the signal expectations (rate, two-dimensional distribution) 
in this channel are obtained by summing the contributions from the three original signals 
weighted by their expected rates. 



7 Higgs boson searches in events with missing energy 
and jets 

The signal topology in this channel is characterized by two acollinear jets and a large 
imbalance in the energy collected by the detector compared to the collision energy, due 
to neutrinos coming either from the decay of a Z boson or from the fusion process. In 
addition to the irreducible qq vv four-fermion background, several other backgrounds can 
lead to similar topologies, like beam-related backgrounds or the qq(7) process with initial 
state radiation photons emitted along the beam axis. Thus, a correct description of the 
initial state radiation in the qq(7) generator particularly matters in this channel. As 
differences of a few % were observed in the a/s 7 distributions obtained with the current 



qq(7) generator [Q and with an analytical calculation [p4]l , the simulated qq(7) events 
have been reweightedQto reproduce the analytical result. This correction has been applied 
throughout the analysis which is described below. 

Events due to particles of the beam with momenta far from the nominal values are 
first excluded by requiring at least two charged particles with impact parameters less 
than 1 mm in the transverse plane and less than 3 mm along the beam direction, and 
with a transverse momentum greater than 2 GeV/c. A loose hadronic preselection is 
then applied, requiring at least nine charged particles, a total charged energy greater 
than 0.16-y/s, a transverse energy greater than 0.15y/s and the sum of the magnitudes of 
all particle momenta resolved along the thrust axis to be greater than 0.25a/s- Finally, 
events with an electromagnetic shower exceeding 0.45a/s are rejected. These criteria 
remove 97% of the 77 background and veto completely the Bhabha background. 

In order to reject events coming from a radiative return to the Z with photons emitted 



in the beam pipe, V V is required to be greater than 115 GeV when the polar angle of 

2 Events with sis' above (resp. below) 0.85-^/s are reweighthed by a factor 0.96 (resp. 1.025). 
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the missing momentum is within 40° to the beam axis. To reduce the contamination of 
radiative return events with photons in the detector acceptance, events are rejected if their 
total electromagnetic energy within 30° to the beam axis is greater than 0.16y/s or if the 
total energy in the small angle luminosity monitor is greater than 0.08-y/s. A veto based 
on the hermeticity counters of DELPHI as described in [ 20 1 is also applied to reject events 
with photons crossing the small insensitive regions of the electromagnetic calorimeters. 
To reduce the two-fermion background outside the radiative return peak as well as four- 
fermion backgrounds without missing energy, y/s* must not exceed 0.96a/s. Two-fermion 
events with jets pointing to the insensitive regions of the electromagnetic calorimeters 
are also a potential background due to mismeasurements of the jet properties. To reject 
such a background, events are forced into a two-jet configuration using the DURHAM 
algorithm and are rejected if the jet polar angles are within ±5° of 40° for one jet and 
of 140° for the other jet, unless the acoplanarity^ is greater than 10°. At this stage, 88% 
of the total qq(7) background is removed. In order to reduce most of the contamination 
from semi-leptonic decays of W + W~ pairs, the energy of the most energetic particle of 
the event must not exceed 0.2^. To reinforce the rejection of those decays containing 
a t lepton, there must be no charged particle in the event with a transverse momentum 
with respect to its jet axis greater than 10 GeV/c when forcing the event into the two-jet 
configuration. The final selection of signal-like events requires the total visible energy to 
be lower than 0.70^. All the above criteria define the preselection. 

The final discrimination between signal and background is achieved through a multi- 
dimensional variable built with the likelihood ratio method. A short description of the 
algorithms needed in this step is given below. As already mentioned, events are forced 
into two jets with the DURHAM algorithm (the so called "two-jet configuration") but 
for each event jets are also reconstructed with the same algorithm using a distance of 
Ucut = 0.005 (the so called "free-jet configuration") and general variables of each jet (like 
multiplicities, momenta) are calculated in both configurations. In order to tag remain- 
ing isolated particles from semi-leptonic decays of W + W~ pairs, the energies collected 
between two cones with half opening angles of 5° and 25° around the most isolated and 
the most energetic particles are calculated and normalised to the corresponding particle 
energies. The lesser of these two normalised energies defines the anti-W + W~ isolation 
variable. 

The likelihood multidimensional variable combines the following discriminant vari- 
ables: the angle between the missing momentum and the closest jet in the free-jet con- 
figuration, the polar angle of the more forward jet in the two-jet configuration, the polar 
angle of the missing momentum, the acoplanarity in the two-jet configuration, the ratio 
between y/s' and the centre-of-mass energy, the missing mass of the event, the anti-W + W~ 
isolation variable, the largest transverse momentum with respect to its jet axis of any 
charged particle in the two-jet configuration, the DURHAM distance for the transition 
betwen the two-jet and three-jet configurations, the minimum jet charged multiplicity in 
the free-jet configuration, the event lifetime probability P£ that all tracks with a posi- 
tive lifetime-signed impact parameter in the event give a product of track significances as 
large as that observed if they do come from the interaction point, and the global b-tagging 
variable x\>- The first five variables discriminate the signal from the qq(7) background 
and the other variables provide a discrimination against W + W~ pairs. For each variable, 
probability density functions (p.d.f.s) at each centre-of-mass energy are obtained from 
simulated events, using half of the statistics available in all backgrounds and in signals of 
masses 95, 100 and 105 GeV/c 2 at below 198 GeV, and 100, 105, 110 and 115 GeV/c 2 

3 The acoplanarity is defined as the supplement of the angle between the transverse momenta (with respect to the beam 
axis) of the two jets. 
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for above 198 GeV. The whole samples are used to derive the final results, in order to 
improve limited statistics in some bins of the two-dimensional discriminant information 
used to derive the confidence levels. 

The distributions of four of the input variables are shown at preselection level in Fig. |5], 
while Fig. |6] shows the distribution of the likelihood discriminant variable. The compari- 
son between the observed and expected rates in the signal-like tail of this distribution is 
illustrated further in Fig. ^ which shows the observed and expected background rates at 
= 199.6 GeV/c 2 as a function of the efficiency for a Higgs signal of 105 GeV/c 2 when 
varying the cut on the likelihood variable. As a final selection, a minimal value of 1.0 is 
required, leaving 108 events in data for a total expected background of 105.7 ± 1.2(stat.). 
The two-dimensional calculation of the confidence levels uses the likelihood variable and 
the recontructed Higgs boson mass defined as the visible mass given by a one-constraint 
fit where the recoil system is an on-shell Z boson. The effect of the selections on data 
and simulated samples are detailed in Tables p] to |4], while the efficiencies at the end of 
the analysis are reported in Tables |5| and ^| as a function of wh ■ 

Systematic uncertainties due to the use of non-independent samples in the definition 
of the p.d.f.s and in the final result derivation were estimated by comparing the results 
when running the analysis on the reference samples and on the complementary samples. 
The differences between the two sets of results are then quoted as systematics if they 
are higher than the statistical uncertainties. These systematics amount to ±2.0% for the 
efficiencies and to ±4.5% for the background estimates at ^/s = 191.6 and 201.7 GeV/c 2 
while no significant difference is observed for the background estimates at the other two 
energies. Systematic uncertainties due to the imperfect modelling of the detector re- 
sponse were derived by rescaling the bin contents of each p.d.f. from simulation to those 
in data, restricting to bins where the deviation between data and simulation exceeded 
two standard deviations. The analysis was then repeated with the rescaled p.d.f. for 
each variable in turn and the largest difference with respect to the initial result taken 
as systematics. These amount to ±2.0% on the efficiencies and ±10.0% on the back- 
ground estimates and come from the p.d.f. of the acoplanarity. It was checked that these 
differences remained similar with tighter selections in the likelihood variable. Thus the 
overall uncertainties are ±3.0% on the efficiencies, ±10.0% (±11.0%) on the background 
estimates at = 195.6 and 199.6 GeV/c 2 (191.6 and 201.7 GeV/c 2 ). 

A second analysis using the same preselection criteria followed by an Iterated nonlinear 
Discriminant Analysis (IDA) as described in |4j] gave similar results. 

8 Higgs boson searches in pure hadronic events 

Higgs boson searches in pure hadronic final states start with a common four-jet pre- 
selection, which eliminates 77 events and reduces the qq(7) and Z7* backgrounds. As 
this step did not change since the previous analysis, the reader is referred to 0,|T|] for 
the exact description of the cuts and only the important features are briefly mentioned 
here. After a selection of multi-hadron events excluding those with an energetic photon 
in the calorimeters or lost in the beam pipe, topological criteria are applied to select 
multi-jet events. All selected events are then forced into a four-jet topology with the 
DURHAM algorithm and a minimal multiplicity and mass is required for each jet. After 
the preselection, different analysis procedures are applied in the HZ and hA channels. 
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8.1 The HZ four-jet channel 

After the common four-jet preselection, events are selected using a discriminant vari- 



able defined as the output of an artificial neural network |25| which combines four vari- 
ables. Three of them are introduced to reduce the four-fermion contamination. The first 
relies on b-tagging and is the maximum b-tagging variable of any di-jet in the event, a 
di-jet b-tagging variable being defined as the sum of the two jet b-tagging variables, x l h . 
The second and third variables rely on mass information and test the compatibility of 
the event with the hypotheses of W + W~ and ZZ pair-production, respectively. First, 
constrained fits are used to derive the probability density function measuring the com- 
patibility of the event kinematics with the production of two objects of any masses. This 
two-dimensional probability, called the ideogram probability, is then folded with the ex- 
pected mass distributions for the W + W~ and ZZ processes, respectively. More about the 
ideogram technique can be found in [p6[| . Finally, the fourth input variable to the neural 
network is intended to reduce the qq(7) contamination and is the output of another neu- 
ral network [^] (anti-QCD neural network) constructed from eight variables. These are 
mostly shape or jet variables: the sum of the second and fourth Fox- Wolfram moments, 
the product of the minimum jet energy and the minimum opening angle between any 
two jets, the maximum and minimum jet momenta, the sum of the multiplicities of the 
two jets with lowest multiplicity and the sum of the masses of the two jets with lowest 
masses. For the last two variables, the six possible pairings of the jets are considered and 
the variables are defined as the minimum di-jet mass and the minimum sum of the cosines 
of the opening angles of the two dijets in any pairing. As the discrimination between the 
qq(7) background and the Hqq signal provided by these variables depends mainly on 
the difference y/s— m H , the anti-QCD neural network was trained with simulations at 
y/s = 189 GeV, using qq(7) events and 95 GeV/c 2 Hqq events. In order to minimize 
the risk of overtraining, the neural network used for the final discrimination was trained 
with fractions of the available simulated samples at y/s = 195.6 GeV in qq(7) background 
(10%), four-fermion background (50%), and 105 GeV/c 2 Hqq signal (50%). The whole 
samples were used to derive the final results. 

The agreement between data and background simulation after the four-jet preselection 
is illustrated in Fig. [7] which shows the distributions of three analysis variables and of the 
recontructed Higgs boson mass obtained as explained below. Fig. [| shows the distribution 
of the final neural network output variable and, as an example, the expected background 
rate and the data at y/s = 199.6 GeV, as a function of the efficiency for a 105 GeV/c 2 
signal when varying the cut on the neural network output variable. As a final selection, a 
minimal value of 0.3 is required. This suppresses the most background-like events, leaving 
161 events in data and a total expected background of 175.4 ± 1.3(stat.). The effect of the 
selections on data and simulated samples are detailed in Tables [I] to || while the efficiencies 
at the end of the analysis are reported in Tables and ^| as a function of mn for the 
Hqq channel and in Table [7] as a function of wh and tua for the (h — > AA) qq channels. 
Since these two channels, specific to the MSSM, are covered by the same analysis as 
that of the Hqq channel, the three channels cannot be considered as independent in the 
confidence level computation when testing MSSM models. For this computation, they 
are combined into one global Hqq channel: at each test point, the signal expectations 
(rate and two-dimensional distribution as defined below) in this channel are obtained by 
summing the contributions from the three original signals weighted by their expected 
rates. 

The two-dimensional calculation of the confidence levels uses the final neural network 
variable and the recontructed H boson mass estimated as follows. For each of the six 



11 



possible pairings of jets into an HZ pair, a kinematic fit is applied, requiring energy and 
momentum conservation and one di-jet to be at the nominal Z mass. The pairing of jets 
defining the Higgs boson and Z candidates is then that which maximises the probability ||] 
that both the b-content of the different jets and the \ 2 probability of the five-constraint 
fit are compatible with the production of an HZ pair. 

The systematic uncertainties from the imperfect modelling of the detector response 
were estimated by repeating the selection procedure on the distribution of the neural 
network variable obtained by smearing, in turn, each of the distributions of the three input 
variables according to the resolution in the variable. This leads to relative uncertainties of 
±6.0% related to b-tagging, ±3.0% related to the anti-QCD variable and ±2.5% related 
to the WW ideogram probability. This results in an overall relative uncertainty of ±7.2% 
in the background and efficiency estimates at each centre-of-mass energy. 

8.2 The hA four-b channel 

The analysis is very similar to that published in After the common four-jet pres- 
election, events are preselected further, requiring a visible energy greater than 120 GeV, 
^fs' greater than 150 GeV, a missing momentum component along the beam direction 
lower than 30 GeV/c and at least two charged particles per jet. A four-constraint kine- 
matic fit requiring energy and momentum conservation is then applied, and the two di-jet 
masses are calculated for each of the three different jet pairings. As the possible produc- 
tion of MSSM Higgs bosons through the hA mode dominates at large tan/3 where the 
two bosons are almost degenerate in mass, the pairing defining the Higgs boson candi- 
dates is chosen as that which minimizes the mass difference between the two di-jets. The 
final discrimination between background and signal is then based on a multidimensional 
variable which combines the following eight variables with a likelihood ratio method: the 
event thrust, the second and fourth Fox- Wolfram moments, the difference between the 
Higgs boson candidate masses as given by the kinematic fit, the production angle of the 
Higgs boson candidates, the sum of the four jet b-tagging variables, the minimum di-jet 
b-tagging variable and the number of secondary vertices. For each variable, probability 
density functions (p.d.f.s) were obtained from simulated events, using fractions of the 
statistics available in the qq(7) background (40%) and four-fermion background (80%) 
at y/s = 195.6 GeV and 199.6 GeV and in signal events with = 85, 90 GeV/c 2 and 
tan/3 =20 (50%) at y/s = 195.6 GeV. The whole samples were used to derive the final 
results. 

The agreement between data and background simulation after the preselection is illus- 
trated in Fig. ^ which shows the distributions of three input variables and of the sum of 
the recontructed Higgs boson masses as given by the kinematic fit. Fig. |TD| shows the dis- 
tribution of the final discriminant variable and, as an example, the expected background 
rate and the data at i/i = 199.6 GeV, as a function of the efficiency for a signal with 
mp = 85 GeV/c 2 and tan/? = 20, when varying the cut on the discriminant variable. As 
a final selection, a minimal value of 0.1 is required, leading to 136 events in data, for a 
total expected background of 137.8 ± 1.2(stat.). The effect of the selections on data and 
simulated samples are detailed in Tables [l] to f|, while the efficiencies at the end of the 
analysis are reported in Tables ||| and ^ as functions of m\ and tan f3 and of wa and m^. 

The two-dimensional calculation of the confidence levels uses the likelihood variable 
and the sum of the reconstructed Higgs boson masses as given by the kinematic fit. 

Systematic uncertainties due to the use of non-independent samples in the definition of 
the p.d.f.s and in the final result derivation were estimated at the level of ±4.0% relative, 
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by repeating the whole procedure with two independent samples of lower size. Systematic 
uncertainties due to the imperfect modelling of the detector response were derived as in 
the previous section. The uncertainty related to b-tagging amounts to ±5.0% and that 
related to shape variables to ±3.0%, resulting in an overall relative uncertainty of ±7.0% 
on background and efficiency estimates at each centre-of-mass energy. 



9 Results 

The results of the searches presented in the previous sections can be translated into 
exclusion limits on the masses of the neutral Higgs bosons in the SM and MSSM. 



9.1 Reconstructed mass spectra 

As an illustration of the discrimination achieved against the residual SM background, 
distributions of the reconstructed Higgs boson mass(es) after tight selections are presented 
in Fig. [II] in the HZ and hA channels. The selections correspond to requiring a minimal 
b-tagging value of -1.8 in the He + e~ and H/i + yU~ channels, minimal likelihood values of 
0.8, 7.0 and 3.5 in the r + r~qq , B.uu and 4b channels, respectively, and a minimal neural 
network output of 0.85 in the Hqq channel. The corresponding observed and expected 



rates at each of the four centre-of-mass energies are summarized in Table 10 



9.2 The SM Higgs boson 

We proceed to set a limit on the SM Higgs boson mass, combining the data analysed 
in the previous sections with those taken at lower energies, namely 161.0, 172.0 GeV ppf , 
182.7 GeV [Q and 188.7 GeV The expected cross-sections and branching ratios are 
taken from the database provided by the LEP Higgs working group, using the HZHA [IB] 
package, Version 3, with the top mass set to 174.3 GeV/c 2 . 

Curves of the confidence level CLb and CL S as a function of the test mass mn are 
shown in Fig. |12[ In the presence of a sizeable Higgs signal, the value of the observed CLb 



(top of Fig. [12]) would approach one, since it measures the fraction of background-only 
experiments which are more background-like than the observation. Here the compati- 
bility between the observation and the expectation from background-only is well within 
one standard deviation over the range of masses tested. Moreover, the mass giving an 
expected 5a discovery, defined by the intersection of the curve for signal plus background 
experiments with the horizontal line at 1 — CLb = 5.7 x 10~ 7 , is 98.2 GeV/c 2 . The 
pseudo-confidence level in the signal is shown in Fig. [12] (bottom). The observed 95% CL 
lower limit on the mass is 107.3 GeV/c 2 while the expected median limit is 106.4 GeV/c 2 . 
The curve of the test-statistic Q as a function of the mass hypothesis is shown in 



Fig. |T3 , where the observation is compared with the expectations from background-only 
experiments (top) and from signal plus background experiments (bottom). Over the 
whole range of masses, the test-statistic remains positive, while in the event of a discovery 
it would be negative for mass hypotheses close to the actual mass of the signal. 



9.3 Cross-section limit 

In a more general approach, the results of the searches for a SM Higgs boson can be 
used to set a 95% CL upper bound on the Higgs boson production cross-section, assuming 
that the Higgs boson decay properties are identical to those in the SM but that the Higgs 
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boson couplings to pairs of Z and W ± bosons (the latter arising in the W + W~ fusion 
production mechanism) may be smaller. To achieve the best sensitivity over the widest 
range of mass hypotheses, the results described in this paper are combined with those 
obtained at lower energies at LEP2 PPPOfl, as well as with those obtained at LEP1 [28 



which covered masses up to 60 GeV/c 2 . Both sets of results are treated with the same 
statistical procedure as for the SM. For each mass hypothesis, the production cross- 
section is decreased with respect to its SM value until a pseudo-confidence level CL S of 
5% is obtained. The result is shown in Fig. [TJ as an upper bound on the production 
cross-section, normalised to that in the SM, for masses of the Higgs boson from to 
110 GeV/c 2 . The SM result described in the previous section corresponds to a ratio of 
1. 



9.4 Neutral Higgs bosons in the MSSM 

The results in the hZ and hA channels reported in the previous sections are com- 
bined with the same statistical method as for the SM, also using earlier results at LEP2 
energies [|,ggU|,gU|. The exclusion limits obtained at LEP1 (m h >44 (46) GeV/c 2 
when mj] is above (below) the AA threshold) are used as external constraints to limit the 
number of points in the scans. 

9.4.1 The benchmark scenarios 

At tree level, the production cross-sections and the Higgs branching fractions in the 
MSSM depend on two free parameters, tan (3 and one Higgs boson mass, or, alternatively, 
two Higgs boson masses, eg and m h . Radiative corrections introduce additional pa- 
rameters, related to supersymmetry breaking. Hereafter, we make the usual assumption 
that some of them are identical at a given energy scale: hence, the SU(2) and U(l) gaug- 
ino mass terms are assumed to be unified at the so-called GUT scale, while the sfermion 
mass terms or the squark trilinear couplings are assumed to be unified at the EW scale. 
Within these assumptions, the parameters beyond tree level are: the top quark mass, the 
Higgs mixing parameter, /i, the common sfermion mass term at the EW scale, M susy , the 
SU(2) gaugino mass term at the EW scale, M 2 , the gluino mass, m g , and the common 
squark trilinear coupling at the EW scale, A. The U(l) gaugino mass term at the EW 
scale, Mi, is related to M 2 through the GUT relation Mi = (5/3)tan 2 6 W M 2 . The radia- 
tive corrections affect the relationships between the masses of the Higgs bosons, with the 
largest contributions arising from the top/stop loops. As an example, the h boson mass, 
which is below that of the Z boson at tree level, increases by a few tens of GeV/c 2 in 
some regions of the MSSM parameter space due to radiative corrections. 



In the following, we consider three benchmark scenarios, as suggested in [31]. The first 



two schemes, called the m™ ax scenario and the no mixing scenario, rely on radiative cor- 



rections computed at two-loop order as in p2j. The values of the underlying parameters 



are quoted in Table [11]. The two scenarios differ only by the value of X t = A — fi cot (3, 
the parameter which controls the mixing in the stop sector, and hence has the largest 
impact on the mass of the h boson. The m™ ax scenario leads to the maximum possible 
h mass as a function of tan f3. The no mixing scenario is its counterpart with vanishing 
mixing, leading to upper bounds on which are at least 15 GeV/c 2 lower than in the 

m max sc ] leme _ 

The third scenario, called the large /i scenario, predicts at least one scalar Higgs boson 
with a mass within kinematic reach at LEP2 in each point of the MSSM parameter space. 
However, there are regions for which the Higgs bosons fall below detectability because 



14 



of vanishing branching fractions into b quarks due to large radiative corrections. In this 



scenario, the radiative corrections are computed as in [33]. The values of the underlying 



parameters are given in Table |Tl|. The main difference with the two previous schemes is 
the large and positive value of fx and the relatively small value of trig. 

9.4.2 The procedure 

In the three benchmark scenarios, a scan is made over the MSSM parameters tan/? 
and 771a- The range in m\ spans from 12 GeV/c 2 , the minimal value which has been 
searched for at LEP2 in the DELPHI analyses, up to the maximal value allowed by each 
scenario |[31|| , that is up to M susy , which is 1 TeV/c 2 in the m™ ax and no mixing schemes, 
and 400 GeV/c 2 in the large fx scenario (see Table [11]). The range in tan/3 goes from the 
minimal value allowed in each scenario (0.7 in the large fx scenario and 0.4 in the other 
two schemes) up to 50, a value chosen in the vicinity of the ratio of the top- and b-quark 
masses, which is an example of the large tan/? hypothesis favored in some constrained 
MSSM models ||34|| . The scan steps are 1 GeV/c 2 in mx and 0.1 in tan/? in the regions 
where varies rapidly with these parameters. 

At each point of the parameter space, the hZ and hA cross-sections and the Higgs 
branching fractions are taken from theoretical databases provided by the LEP Higgs 



working group |J5] on the basis of the theoretical calculations in |32|j33|1 . The signal 
expectations in each channel are then derived from the theoretical cross-sections and 
branching fractions, the experimental luminosity and the efficiencies. A correction is 
applied to account for differing branching fractions of the Higgs bosons into bb and t + t~ 
between the test point and the simulation (e.g. for the hZ process, the simulation is done 
in the SM framework). For the hA channels, to account for the difference between the 
masses of the h and A bosons at low tan j3 as well as for the non-negligible width of the 
h and A bosons at large tan/3, the set of efficiencies as a function of obtained from 
the simulations at tan /? = 50 are applied above 30 in tan /3, while the efficiencies derived 
from the tan/? = 20 (or tan/? = 2) simulations are applied between 2.5 and 30 (or below 
2.5) provided the difference between and wa at the test point is below 25 GeV/c 2 ; 
otherwise the set of efficiencies as a function of and tua derived from the additional 
simulations corresponding to large mass differences between the two bosons is preferred. 
The same holds for the discriminant information. Finally, as there is a large overlap in 
the backgrounds selected by the analyses in the Hqq and 4b channels, only one channel 
is selected at each input point and at each centre-of-mass energy, on the basis of the best 
expected CL S from background-only experiments. This ensures that the channels which 
are then combined in the global confidence level computations are independent. 

9.4.3 Results 

To illustrate the compatibility tests of data with background only and with signal plus 
background hypotheses in the hA channels, Fig. ^ shows the curves of the test-statistic 
Q and of the confidence levels CLb and CL S as a function of the test mass mi,+mA, 
when using only the results in the two hA channels. The signal cross-sections are from 
the m™ ax scenario at tan/? around 20. Over the whole range of test masses, data are in 
reasonable agreement with the background expectations. The largest deviation, slightly 
over one standard deviation, is observed for test masses m h +mA around 135 GeV/c 2 and 
is due to the small excess of events in the 4b channel with reconstructed masses in that 



region, as seen in Fig. |TT . 
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Combining the results in the hZ and hA channels gives regions of the MSSM parameter 
space which are excluded at 95% CL or more. The excluded regions in the (mh, tan/3), 
(mA, tan/3) and (mh, mx) planes are presented in Fig. [16] for the m™ ax scenario and in 



Fig. |T7|for the no mixing scenario. Basically, the exclusion is made by the results in the hZ 
(hA) channels in the low (large) tan/3 region while they both contribute at intermediate 
values. For below the kinematic threshold = 2mA, which occurs at low tan (3 only, 
the decay h— >AA opens, in which case it supplants the h— >bb decay. However, in most 
of the region, the A— >bb branching fraction remains large which explains why the results 
in the two (h — > AA) qq channels reported in section |S.1| , combined with studies of the 
h — ► AA decay at lower energies [fHj20|, exclude most of this region. An unexcluded hole 
remains in the no mixing scenario at tan/3 ~ 0.4, between 20 and 40 GeV/c 2 and m^ 
around 85 GeV/c 2 (visible only in the (m\, tan/3) and (mh, wa) projections). In that 
area, the A— >cc decay dominates over the A— >bb decay but the branching fractions in 
both modes are no longer large enough to give the necessary sensitivity for an exclusion. 

The above results establish 95% CL lower limits on and mA, for either assumption 
on the mixing in the stop sector and for all values of tan /3 above 0.49: 

m h > 85.9 GeV/c 2 m A > 86.5 GeV/c 2 . 

The expected median limits are 86.4 GeV/c 2 for and 87.0 GeV/c 2 for m^. The limit 
in mA is reached in the no mixing scenario at tan/? around 30 and thus is due to the 
non-negligible widths of the Higgs bosons, while the limit in is obtained in the m™ ax 
scenario at tan /3 around 7, in a region where both the hZ and hA processes contribute. 
Furthermore, there are excluded ranges in tan/3 between 0.49 and 3.86 (expected [0.49- 
3.86]) in the no mixing case and between 0.65 and 1.75 (expected [0.72-1.75]) in the m™ ax 
scenario. 

The excluded regions in the large /1 scenario are presented in the (mh, tan/3) and (mA, 



tan 0) planes in Fig. [18]. A large fraction of the allowed domain is excluded by the present 
results in the hZ and hA channels. In particular, given that the theoretical upper bound 
on the h boson mass in that scenario is slightly above 107 GeV/c 2 , the sensitivity of the 
hZ channels is high even at large tan /3, which explains why the excluded region reaches 
the theoretically forbidden area for values of tan/3 up to 13.5. On the other hand, there 
is an unexcluded hole in the low tan/3 region at mh around 60 GeV/c 2 which is due to a 
loss of sensitivity because of vanishing h— >bb branching fractions in that region. 



9.4.4 Extended scan of the parameter space 

The robustness of the limits obtained in the benchmark scenarios has been tested in an 
extended scan of the MSSM parameter space. The Higgs bosons masses, cross-sections 
and branching fractions are computed with radiative corrections at two-loop order as 
in [[!2J]. The top mass is fixed at 175 GeV/c 2 while the MSSM parameters, mA, tan/3 
and the parameters governing the radiative corrections, M susy , M 2 , /x and A are varied 
within the ranges given in Table [12]. The values /i = ±1000 GeV/c 2 have been studied 
in addition. As far as the granularity of the scan is concerned, steps of 1 GeV/c 2 are 
used for uia up to 200 GeV/c 2 and larger steps between 200 and 1000 GeV/c 2 ; for each 
value of mA, up to 2700 parameter combinations are investigated. To limit the number of 
points in the scan, only points above 70 GeV/c 2 in and mA are considered, since all 
points below this limit have already been excluded by our previous extended search |36| . 
The scan relies on the same channels and data sets as the representative scans previously 
reported and uses the same procedure to compute the confidence levels at each input 
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point. However, for some parameter sets, the branching ratio of the neutral Higgs bosons 
into neutralinos is dominant, which is never the case in the benchmark scenarios. In such 
a case, the cross-section limits obtained in the search for invisible decays of a neutral 
Higgs boson [37] are applied to check whether these points are excluded or not. Any 
other point in a given plane (e.g. the (mi, tua) plane) is excluded if the observed CL S 
at that point is below 5% for all sets of values of the parameters governing the radiative 
corrections that correspond to that point. The results of this extended scan are presented 
in Fig. O in the three projections (m^, tan/3), (tua, tan/3) and (m^, m^). The extension 
of the MSSM parameter ranges in the scan leads to 95% CL lower limits of 85 GeV/c 2 
on rrih and 86 GeV/c 2 on tua, thus only about 1 GeV/c 2 below the limits obtained in the 
m max an( j no m j x j n g scenarios. 



10 Conclusions 

The 228 pb" 1 of data taken by DELPHI at 191.6-201.7 GeV, combined with our lower 
energy data, sets the lower limit at 95% CL on the mass of the Standard Model Higgs 
boson at: 

m H > 107.3 GeV/c 2 . 

These data sets also allow studies of the representative m™ ax and no mixing scenarios. 
The 95% CL limits on the masses of the lightest neutral scalar and neutral pseudoscalar 
are: 

m h > 85.9 GeV/c 2 m A > 86.5 GeV/c 2 . 

for all values of tan/3 above 0.49 and assuming rriA> 12 GeV/c 2 . These limits have been 
proved to be robust in an extended scan of the MSSM parameter space. 
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Selection 


Data Total qq(7) 4 fermion Efficiency (%) 




background 




Electron channel 25.2 pb 1 


Preselection 


152 158.6± 2.5 111.7 43.6 


81.0 


cuts on leptons 


15 9.7± 0.6 2.5 6.2 


65.7 


5C fit prob. 


3 3.7± 0.3 1.1 2.3 


61.7 


final selection 


1 1.19± 0.08 0.05 1.1 


55.0 


Muon channel 25.9 pb 1 


Preselection 


336 364.8 ± 3.8 274.5 86.7 


78.8 


cuts on leptons 


2 2.30 ± 0.12 0.15 2.15 


72.4 


final selection 


1 0.93 ± 0.02 0.0 0.93 


57.8 


Tau channel 25.9 pb 1 


Preselection 


1209 1127± 2.6 747 380 


98.3 




2.1± 0.1 0.09 2.0 


18.5 


final selection 


0.76±0.08 0.04 0.72 


16.7 


Missing energy channel 24.9 pb 1 


Anti 77 


2378 2368.0± 4.5 1904.1 427.3 


86.0 


Preselection 


139 130.7 ± 2.0 81.3 45.7 


73.9 


C> 1.0 


11 12.6± 0.6 7.7 4.9 


59.8 


Four-jet channel 25.9 pb 1 


Preselection 


302 280.6± 2.8 91.8 188.8 


89.4 


ANN > 0.3 


16 19. 1± 0.5 4.3 14.8 


60.4 


hA four-jet channel 25.9 pb 1 


Preselection 


273 255.1± 2.3 79.7 175.4 


90.7 


C> 0.1 


18 16.8 ± 0.6 7.7 9.1 


85.0 



Table 1: Effect of the selection cuts on data, simulated background and simulated signal 
events at y/s = 191.6 GeV. Efficiencies are given for a signal with m H = 105 GeV/c 2 for 
the SM and m A = 85 GeV/c 2 , tan/3 = 20 for the MSSM. The quoted errors are statistical 
only. For each channel, the first line shows the integrated luminosity used; the last line 
gives the inputs for the limit derivation. 
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Selection 


Data Total qq(7) 4 fermion Efficiency (%) 




background 




Electron channel 76.2 pb 1 


Preselection 


452 428.8± 3.5 284.7 133.7 


79.1 


cuts on leptons 


33 29.0± 1.3 7.7 18.5 


64.2 


5C fit prob. 


18 11.5± 0.7 3.15 7.5 


60.6 


final selection 


5 3.88±0.18 0.13 3.5 


57.4 


Muon channel 76.9 pb 1 


Preselection 


1081 1092.5 ± 4.2 801.5 280.2 


79.5 


cuts on leptons 


3 7.19 ± 0.17 0.46 6.73 


71.4 


final selection 


2 3.02 ± 0.07 0.02 3.0 


67.4 


Tau channel 76.9 pb 1 


Preselection 


3479 3215 ± 4.7 2056 1159 


98.3 




7 6.6± 0.2 0.5 6.1 


18.5 


final selection 


3 2.38±0.12 0.07 2.31 


18.3 


Missing energy channel 75.0 pb 1 


Anti 77 


7005 6757.8 ± 5.2 5343.0 1309.1 


86.1 


Preselection 


403 384.7 ± 2.2 238.8 134.1 


74.4 


£> 1.0 


38 34.0 ± 0.6 20.2 13.8 


62.1 


Four- jet channel 76.9 pb 1 


Preselection 


839 827.1± 5.1 260.0 567.1 


87.5 


ANN > 0.3 


51 58. 7± 0.7 11.3 47.1 


67.7 


hA four-jet channel 76.9 pb 1 


Preselection 


747 757.3± 4.7 232.3 525.0 


90.9 


£> 0.1 


47 48.0± 0.7 19.6 28.5 


86.8 



Table 2: As in Table [I], but for ^ = 195.6 GeV. 
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Selection 


Data Total qq(7) 4 fermion Efficiency (%) 




background 




Electron channel 82.8 pb 1 


Preselection 


489 453.3± 3.6 294.5 148.4 


79.3 


cuts on leptons 


30 31. 1± 1.4 7.0 21.3 


63.3 


5C fit prob. 


11 12.8±1.1 2.6 8.4 


59.5 


final selection 


4 4.22±0.19 0.05 3.9 


56.0 


Muon channel 84.3 pb 1 


Preselection 


1141 1148.4 ± 4.3 807.3 329.8 


78.9 


cuts on leptons 


11 8.20 ± 0.21 0.54 7.66 


72.1 


final selection 


5 3.59 ± 0.08 0.02 3.57 


69.7 


Tau channel 84.3 pb 1 


Preselection 


3629 3434 ± 9.7 2152 1282 


98.1 




8 7.7± 0.2 0.5 7.2 


18.9 


final selection 


3 2.60±0.13 0.14 2.46 


17.9 


Missing energy channel 82.2 pb 1 


Anti 77 


7211 7112.4± 5.8 5566.8 1450.8 


85.5 


Preselection 


421 425.7 ± 2.5 260.6 151.8 


75.8 


£> 1.0 


38 40.5± 0.7 24.4 16.1 


64.0 


Four- jet channel 84.3 pb 1 


Preselection 


882 896.8± 2.9 273.6 623.3 


87.7 


ANN > 0.3 


61 65.3± 0.8 13.1 52.2 


70.4 


hA four-jet channel 84.3 pb 1 


Preselection 


783 817.9± 2.7 243.1 574.8 


90.9 


£> 0.1 


44 49. 1± 0.7 18.9 30.2 


85.5 



Table 3: As in Table [I], but for ^ = 199.6 GeV. 



22 



Selection 


Data Total qq(7) 4 fermion Efficiency (%) 




background 




Electron channel 40.4 pb 1 


Preselection 


232 214.3± 2.0 136.5 72.7 


79.6 


cuts on leptons 


18 15. 7± 0.7 3.83 10.6 


64.3 


5C fit prob. 


3 6.47± 0.55 1.51 4.1 


60.2 


final selection 


1 2.18±0.10 0.09 1.97 


56.8 


Muon channel 41.1 pb 1 


Preselection 


574 561.9 ± 2.9 391.1 165.2 


81.0 


cuts on leptons 


4.31 ± 0.14 0.28 4.03 


73.5 


final selection 


1.83 ± 0.04 0.0 1.83 


71.0 


Tau channel 41.1 pb 1 


Preselection 


1716 1648± 5.8 1019 629 


98.5 




3.6± 0.1 0.2 3.4 


23.1 


final selection 


1.17± 0.05 0.03 1.14 


22.1 


Missing energy channel 40.4 pb 1 


Anti 77 


3305 3401.1 ± 3.7 2632.3 715.3 


85.8 


Preselection 


209 204.9 ± 1.7 125.2 73.3 


77.8 


£> 1.0 


21 18.6 ± 0.5 11.0 7.5 


65.7 


Four- jet channel 41.1 pb 1 


Preselection 


442 432.4± 3.2 129.8 302.6 


87.0 


ANN > 0.3 


33 32.3± 0.5 6.5 25.8 


69.7 


hA four-jet channel 41.1 pb 1 


Preselection 


405 393.9± 2.1 115.6 278.3 


90.4 


£> 0.1 


27 23.9± 0.4 8.9 14.9 


84.6 



Table 4: As in Table [I], but for ^ = 201.7 GeV. 



23 



m H 
(GeV/c 2 ) 


Electron 
channel 


Muon 
channel 


Ht + t t + t Z 
channel channel 


Mis. Energy 
channel 


Four-jet 
channel 


y/s = 191.6 GeV 



70.0 56.1 ± 1.1 65.0 ± 1.1 21.6 ± 0.9 22.5 ± 0.9 49.0 ± 1.1 45.1 ± 1.1 

80.0 57.4 ± 1.1 70.6 ± 1.0 19.6 ± 0.9 24.0 ± 1.0 57.0 ± 1.1 57.5 ± 1.1 

85.0 59.4 ± 0.8 68.8 ± 1.0 19.6 ± 0.9 23.3 ± 0.9 61.9 ± 1.1 62.3 ± 1.1 

90.0 58.4 ± 0.8 69.0 ± 1.0 20.6 ± 0.9 22.8 ± 0.9 65.0 ± 1.1 69.1 ± 1.0 

95.0 58.1 ± 0.8 70.0 ± 0.8 19.8 ± 0.9 22.0 ± 0.9 66.5 ± 1.1 69.1 ± 1.0 

100.0 55.2 ± 0.8 67.4 ± 1.0 19.5 ± 0.9 19.9 ± 0.9 61.9 ± 1.1 69.4 ± 1.0 

105.0 55.0 ± 0.8 57.8 ± 1.1 16.7 ± 0.8 19.3 ± 0.9 59.8 ± 1.1 60.5 ± 1.1 

110.0 52.5 ± 1.1 47.4 ± 1.1 14.0 ± 0.8 18.4 ± 0.9 59.1 ± 1.1 53.7 ± 1.1 



y/s = 195.6 GeV 



70.0 53.9 ± 1.1 63.3 ± 1.1 20.0 ± 0.9 23.0 ± 0.9 47.9 ± 1.1 42.7 ± 1.1 

80.0 58.3 ± 1.1 68.3 ± 1.0 22.3 ± 0.9 21.8 ± 0.9 55.3 ± 1.1 56.4 ± 1.1 

85.0 59.8 ± 1.1 71.1 ± 1.0 21.3 ± 1.3 24.4 ± 1.4 60.8 ± 1.1 60.4 ± 1.1 

90.0 57.9 ± 1.1 67.7 ± 1.1 20.8 ± 1.3 23.6 ± 1.3 64.2 ± 1.1 66.1 ± 1.0 

95.0 59.8 ± 1.1 69.7 ± 1.1 21.3 ± 1.3 23.4 ± 1.3 67.8 ± 1.0 68.5 ± 1.0 

100.0 59.0 ± 0.9 71.3 ± 0.8 18.2 ± 1.2 21.8 ± 1.3 65.8 ± 0.7 70.0 ± 0.6 

105.0 57.4 ± 1.1 67.4 ± 1.1 18.3 ± 1.2 20.2 ± 1.3 62.1 ± 1.1 67.7 ± 1.0 

110.0 56.0 ± 1.1 55.5 ± 1.1 14.4 ± 0.8 20.4 ± 0.9 55.6 ± 1.1 57.3 ± 1.1 

115.0 53.5 ± 1.1 45.9 ± 1.1 13.8 ± 0.8 16.9 ± 0.8 53.2 ± 1.1 52.4 ± 1.1 



sfs = 199.6 GeV 



70.0 54.6 ± 1.1 62.5 ± 1.1 20.0 ± 0.9 23.3 ± 0.9 44.9 ± 1.1 44.2 ± 1.1 

80.0 56.9 ± 1.1 68.1 ± 1.0 20.7 ± 0.9 23.6 ± 0.9 52.5 ± 1.1 55.1 ± 1.1 

85.0 55.8 ± 1.6 70.2 ± 1.0 22.0 ± 1.5 23.6 ± 1.5 59.5 ± 1.1 59.5 ± 1.1 

90.0 59.6 ± 1.1 69.5 ± 1.0 19.2 ± 1.4 23.9 ± 1.5 62.2 ± 1.1 63.6 ± 1.0 

95.0 59.6 ± 1.1 70.9 ± 1.1 19.9 ± 1.4 21.7 ± 1.5 64.7 ± 1.1 67.3 ± 1.0 

100.0 57.9 ± 1.1 72.1 ± 1.0 20.0 ± 1.4 22.1 ± 1.5 67.7 ± 1.0 69.0 ± 1.0 

105.0 56.0 ± 0.9 69.7 ± 0.6 17.9 ± 1.4 19.5 ± 1.4 64.0 ± 1.0 70.4 ± 0.6 

110.0 56.1 ± 1.1 66.7 ± 1.0 18.3 ± 1.4 20.7 ± 1.4 61.7 ± 0.8 68.2 ± 1.0 

115.0 53.5 ± 1.1 58.3 ± 1.2 14.4 ± 1.1 18.8 ± 1.2 60.4 ± 0.9 60.2 ± 1.5 



y/s = 201.7 GeV 



70.0 53.9 ± 1.1 62.5 ± 1.1 19.6 ± 0.9 24.2 ± 1.0 41.2 ± 1.1 42.5 ± 1.1 

80.0 57.1 ± 1.1 67.2 ± 1.0 19.1 ± 0.9 24.2 ± 1.0 50.5 ± 1.1 53.0 ± 1.1 

85.0 61.4 ± 1.1 69.5 ± 1.1 22.0 ± 1.3 23.6 ± 1.3 59.1 ± 1.1 58.5 ± 1.1 

90.0 58.6 ± 1.1 68.3 ± 1.1 23.0 ± 1.3 23.2 ± 1.3 61.9 ± 1.1 64.5 ± 1.1 

95.0 60.2 ± 1.1 70.0 ± 1.1 21.2 ± 1.3 21.4 ± 1.3 64.3 ± 1.1 63.1 ± 1.1 

100.0 60.2 ± 1.1 67.5 ± 1.1 21.7 ± 1.2 22.0 ± 1.3 66.9 ± 1.0 69.9 ± 1.0 

105.0 56.8 ± 0.9 71.0 ± 1.0 22.1 ± 1.2 22.8 ± 1.3 65.7 ± 1.1 69.7 ± 1.1 

110.0 56.9 ± 1.1 69.1 ± 1.0 15.0 ± 1.1 19.3 ± 1.2 62.1 ± 1.1 69.4 ± 1.0 

115.0 57.3 ± 1.1 58.3 ± 1.2 14.7 ± 1.1 18.8 ± 1.2 63.4 ± 1.1 61.6 ± 1.1 



Table 5: HZ channels: efficiencies (in %) of the selection at \/s = 191.6-201.7 GeV as a 
function of the mass of the Higgs boson, for masses above 70 GeV/c 2 . The quoted errors 
are statistical only. 
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m H (GeV) 
(GeV/c 2 ) 191.6 195.6 199.6 201.7 
Electron channel 

5O0 5.8 ± 0.5 8.5 ± 0.6 8.8 ± 0.6 8.8 ± 0.6 

60.0 37.0 ± 1.1 37.0 ± 1.1 36.2 ± 1.1 38.6 ± 1.1 

Muon channel 

4O0 27.8 ± 1.0 25.0 ± 1.0 23.2 ± 0.9 23.4 ± 0.9 

50.0 44.4 ± 1.1 41.6 ± 1.1 38.4 ± 1.1 37.4 ± 1.1 

60.0 55.0 ± 1.1 55.2 ± 1.1 52.8 ± 1.1 49.4 ± 1.1 

Hr + r" channel 

5O0 12.5 ± 0.7 11.9 ± 0.7 12.7 ± 0.7 12.6 ± 0.7 

60.0 18.9 ± 0.9 18.4 ± 0.9 18.8 ± 0.9 19.1 ± 0.9 

t + t~Z channel 

5O0 7.9 ± 0.6 6.0 ± 0.5 6.3 ± 0.5 6.0 ± 0.5 

60.0 21.2 ± 0.9 17.3 ± 0.8 14.1 ± 0.8 14.1 ± 0.8 

Missing Energy channel 

12T0 17.6 ± 0.8 21.7 ± 0.9 25.2 ± 1.0 18.1 ± 0.9 

18.0 23.8 ± 0.9 32.8 ± 1.0 34.4 ± 1.1 22.4 ± 0.9 

24.0 25.3 ± 0.9 32.9 ± 1.0 33.7 ± 1.1 22.2 ± 0.9 

30.0 25.6 ± 1.0 32.8 ± 1.0 31.5 ± 1.0 23.5 ± 0.9 

40.0 29.0 ± 1.0 32.0 ± 1.0 32.0 ± 1.0 23.3 ± 0.9 

50.0 33.5 ± 1.0 35.7 ± 1.1 34.0 ± 1.1 28.5 ± 1.0 

60.0 36.0 ± 1.1 41.7 ± 1.1 38.1 ± 1.1 35.5 ± 1.1 
Four-jet channel 



24.0 


6.8 ± 0.6 


8.0 ± 0.6 


7.9 ± 0.6 


8.6 ± 0.7 


30.0 


13.5 ± 0.8 


16.2 ± 0.8 


16.6 ± 0.8 


17.1 ± 0.8 


40.0 


29.8 ± 1.0 


29.6 ± 1.0 


28.1 ± 1.0 


28.9 ± 1.0 


50.0 


38.6 ± 1.1 


39.6 ± 1.1 


37.8 ± 1.1 


38.8 ± 1.0 


60.0 


42.9 ± 1.1 


44.0 ± 1.1 


42.4 ± 1.1 


42.1 ± 1.1 



Table 6: HZ channels: efficiencies (in %) of the selection at y/s = 191.6-201.7 GeV as a 
function of the mass of the Higgs boson, for masses below 70 GeV/c 2 . Only efficiencies 
higher than 5% are shown. The quoted errors are statistical only. 
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A^bb A^cc 
wa m h Efficiency Efficiency 



(GeV/c 2 ) 


(GeV/c 2 ) 


(%) 




(%) 




12.0 


30.0 


20.8 


± 


1.2 


6.5 ± 0.8 


12.0 


50.0 


48.5 


± 


1.6 


19.2 ± 


1.2 


12.0 


70.0 


56.0 


± 


1.6 


24.4 ± 


1.4 


12.0 


90.0 


80.2 


± 


1.3 


37.6 ± 


1.5 


12.0 


105.0 


77.3 


± 


1.3 


54.3 ± 


1.6 


20.0 


50.0 


46.0 


± 


1.6 


17.4 ± 


1.2 


20.0 


70.0 


57.6 


± 


1.6 


24.3 ± 


1.3 


20.0 


90.0 


75.6 


± 


1.4 


36.4 ± 


1.5 


20.0 


105.0 


80.1 


± 


1.3 


59.5 ± 


1.5 


30.0 


70.0 


63.4 


± 


1.6 


26.0 ± 


1.4 


30.0 


90.0 


72.8 


± 


1.5 


31.8 ± 


1.5 


30.0 


105.0 


81.4 


± 


1.3 


55.6 ± 


1.6 


40.0 


90.0 


74.3 


± 


1.5 


35.0 ± 


1.6 


40.0 


105.0 


80.3 


± 


1.3 


42.4 ± 


1.6 


50.0 


105.0 


82.4 


± 


1.3 


47.8 ± 


1.7 



Table 7: (h — > AA)(Z — > qq ) channels with A^bb or A^cc : efficiencies of the selection 
(in %) at y/s = 199.6 GeV as a function of the masses of the A and h bosons. The quoted 
errors are statistical only. 





tan/5 


= 2 


tan/5 


= 20 


tan/5 


= 50 


m A 


Four-jet 


Tau 


Four- jet 


Tau 


Four-jet 


Tau 


(GeV/c 2 ) 


channel 


channel 


channel 


channel 


channel 


channel 



y/s = 191.6 GeV 



39.9 ± 1.1 

58.3 ± 1.1 
66.1 ± 1.0 

72.4 ± 1.0 
77.1 ± 0.9 
78.8 ± 0.9 

83.1 ± 0.8 

82.2 ± 0.8 



7.1 ± 
14.6 ± 

22.3 ± 

23.4 ± 
22.6 ± 
20.6 ± 



0.6 
0.8 
0.9 
0.9 
0.9 
0.9 



21.8 ± 0. 
54.8 ± 1. 
67.3 ± 1. 
75.2 ± 1 
81.1 ± 
85.0 ± 0. 
82.8 ± 0. 
78.5 ± 0. 



9 
1 

6.2 ± 0.5 

22.7 ± 0.9 

9 24.3 ± 1.0 

.8 24.1 ± 0.9 
22.1 ± 0.9 

.9 19.0 ± 0.9 



14.5 ± 0.8 

49.5 ± 1.1 

63.6 ± 1.1 
73.4 ± 1.0 
77.0 ± 0.9 



6.5 ± 0.5 
16.9 ± 0.8 
19.9 ± 0.9 



80.8 ± 0.9 20.7 ± 0.9 
79.0 ± 0.9 19.7 ± 0.9 
77.4 ± 1.0 15.4 ± 0.8 



y/i = 195.6 GeV 



40.0 


38.4 ± 1.1 




50.0 


58.0 ± 1.1 




60.0 


67.0 ± 1.0 


7.0 ± 0.6 


70.0 


71.1 ± 1.1 


13.4 ± 0.8 


80.0 


81.0 ± 0.8 


15.9 ± 1.1 


85.0 


81.9 ± 0.8 


22.2 ± 1.3 


90.0 


81.5 ± 0.8 


25.3 ± 1.4 


95.0 


82.6 ± 0.8 


20.6 ± 0.9 



17.7 ± 

55.4 ± 1 
66.9 ± 1 
72.6 ± 1 
83.9 ± 

86.8 ± 
87.0 ± 

82.5 ± 



.1 

.0 7.2 ± 0.6 

.0 21.4 ± 0.9 

8 20.3 ± 1.3 

8 21.8 ± 1.3 

6 22.1 ± 1.3 

8 19.0 ± 0.9 



13.6 ± 0.8 
47.1 ± 1.1 
62.0 ± 1.1 
71.3 ± 1.0 



6.0 ± 0.5 
17.5 ± 0.8 



77.4 ± 0.9 23.7 ± 0.9 

80.9 ± 0.9 21.1 ± 0.9 

82.0 ± 0.9 20.1 ± 0.9 

81.2 ± 0.9 17.6 ± 0.8 



y/s = 199.6 GeV 



38.6 ± 1.1 

58.3 ± 1.1 

67.0 ± 1.0 

71.1 ± 1.1 

78.4 ± 0.9 

80.2 ± 0.8 

83.5 ± 0.8 
84.5 ± 0.7 



5.8 ± 
11.4 ± 
15.9 ± 
20.8 ± 
24.4 ± 
20.6 ± 



0.5 
0.7 
1.1 
1.3 
1.4 
1.3 



14.5 ± 

54.8 ± 1 

66.9 ± 1. 

72.6 ± 1. 
82.2 ± 0. 

85.5 ± 0. 

85.6 ± 0. 
83.5 ± 0. 



1 

6.8 ± 0.6 

19.1 ± 0.9 

9 22.5 ± 1.3 

8 25.6 ± 1.4 

8 21.6 ± 1.3 

8 22.7 ± 1.3 



11.2 ± 0.7 
47.7 ± 1.1 
62.0 ± 1.1 

71.3 ± 1.0 



5.5 ± 0.5 
15.1 ± 0.8 



77.4 ± 0.9 20.1 ± 0.9 

80.9 ± 0.9 20.6 ± 0.9 

82.0 ± 0.9 21.1 ± 0.9 

81.2 ± 0.9 18.9 ± 0.9 



y/s = 201.7 GeV 



33.3 ± 1.0 
60.1 ± 1.1 

67.0 ± 1.0 

71.1 ± 1.1 
77.0 ± 0.9 

80.2 ± 0.8 
82.5 ± 0.8 
85.8 ± 0.8 



5.8 ± 
11.2 ± 
15.9 ± 
20.8 ± 
24.4 ± 
20.6 ± 



0.5 
0.7 
1.1 
1.3 
1.4 
1.3 



13.6 ± 
54.6 ± 1 
66.9 ± 1. 
72.6 ± 1. 
80.4 ± 0. 
84.6 ± 0. 
86.8 ± 0. 
84.1 ± 0. 



1 

6.7 ± 0.6 

17.0 ± 0.8 

9 22.5 ± 1.3 

8 25.6 ± 1.4 

8 21.6 ± 1.3 

8 22.7 ± 1.3 



10.5 ± 0.7 
47.1 ± 1.1 
62.0 ± 1.1 

71.3 ± 1.0 

77.4 ± 0.9 



5.8 ± 0.5 

15.4 ± 0.8 

20.5 ± 0.9 



80.9 ± 0.9 22.4 ± 0.9 
82.0 ± 0.9 20.8 ± 0.9 
81.2 ± 0.9 18.4 ± 0.9 



Table 8: hA channels: efficiencies of the selection (in %) at ^/s = 191.6-201.7 GeV 
a function of the mass of the A boson and tan/5. Only efficiencies higher than 5% ; 
shown. The quoted errors are statistical only. 
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m A 


m h 


Efficiency 


m A 


m h 


Efficiency 


(GeV/c 2 ) (GeV/c 2 ) 


(%) 


(GeV/c 2 ) (GeV/c 2 ) 


(%) 


12.0 


70.0 


27.2 ± 1.4 


30.0 


130.0 


67.7 ± 1.5 


12.0 


90.0 


35.8 ± 1.5 


30.0 


150.0 


60.1 ± 1.5 


12.0 


110.0 


59.8 ± 1.5 


50.0 


70.0 


66.2 ± 1.5 


12.0 


130.0 


63.4 ± 1.5 


50.0 


90.0 


78.4 ± 1.3 


12.0 


150.0 


52.0 ± 1.6 


50.0 


110.0 


80.5 ± 1.3 


12.0 


170.0 


32.4 ± 1.5 


50.0 


130.0 


75.5 ± 1.4 


30.0 


50.0 


19.9 ± 1.3 


70.0 


90.0 


80.4 ± 1.3 


30.0 


70.0 


54.7 ± 1.6 


70.0 


110.0 


81.8 ± 1.2 


30.0 


90.0 


65.2 ± 1.5 


80.0 


90.0 


87.2 ± 1.1 


30.0 


110.0 


68.9 ± 1.5 


85.0 


95.0 


83.5 ± 1.2 



Table 9: hA four-jet channel : efficiencies of the selection (in %) at ^/s = 199.6 GeV as 
a function of the masses of the A and h bosons, from simulated samples corresponding 
to large mass differences between the two bosons. The quoted errors are statistical only. 





HZ channel 




hA channel 


(GeV) Data Background 


Signal 


Data 


Background 


Signal 


191.6 1 


4.8±0.3 


0.16±0.01 


3 


1.1±0.2 


0.56±0.02 


195.6 13 


12.4±0.3 


1.18±0.02 


1 


2.9±0.2 


2.00±0.04 


199.6 13 


12.9±0.3 


4.70±0.06 


3 


2.9±0.2 


2.31±0.05 


201.7 5 


6.9±0.2 


2.87±0.06 


1 


1.5±0.1 


1.24±0.02 



Table 10: Observed and expected rates after tight selections applied to data at 
t/s = 191.6-201.7 GeV. Signal expectations are given for a signal with mn = 105 GeV/c 2 
for the SM and m A = 85 GeV/c 2 , tan/3 = 20 for the MSSM. The quoted errors are sta- 
tistical only. 



scenario 


m top 

(GeV/c 2 ) 


M-susy 

(GeV/c 2 ) 


M 2 
(GeV/c 2 ) 


m- g 
(GeV/c 2 ) 


n 

(GeV/c 2 ) 


X t = A — fj, cot (3 
(GeV/c 2 ) 


m max scenar j 


174.3 


1000 


200 


800 


-200 


2 ^^susy 


no mixing 


174.3 


1000 


200 


800 


-200 





large \x 


174.3 


400 


400 


200 


1000 


-300 



Table 11: Values of the underlying parameters for the three representative MSSM 
scenarios scanned in this paper. 



Parameter 


m A 
(GeV/c 2 ) 


tan f3 


M-susy 

(GeV/c 2 ) 


M 2 
(GeV/c 2 ) 


n 

(GeV/c 2 ) 


Aj M SU gy 


Range 


20: 1000 


0.5:50 


200 : 1000 


200: 1000 


-500: +500 


-2: +2 



Table 12: Ranges of variation of the underlying parameters used in the extended scan 
of the MSSM parameter space described in this paper. 
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Figure 1: b-tagging: Top: distributions of the combined b-tagging variable Xb, in 1999 Z 
data (dots) and simulation (histogram). The contribution of udsc-quarks is shown as the 
dark histogram. Middle: ratio of integrated tagging rates in Z data and simulation as a 
function of the cut in x\>. Bottom: mistag probabilities for c or uds-quarks as a function 
of the efficiency for b-quarks, estimated from simulated Z data. 
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Slow electron energy (GeV) 





Figure 2: He + e~ channel: distributions of four analysis variables, as described in the 
text, at preselection level. Data at yfs = 191.6-201.7 GeV (dots) are compared with SM 
background expectations (left-hand side histograms) and with the expected distribution 
for a 105 GeV/c 2 signal (right-hand side histogram, normalised to 50 times the expected 
rate). 
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Figure 3: H/i + /i~ channel: distributions of four analysis variables, as described in the 
text, at preselection level. Data at y/s = 191.6-201.7 GeV (dots) are compared with SM 
background expectations (left-hand side histograms) and with the expected distribution 
for a 105 GeV/c 2 signal (right-hand side histogram, normalised to one thousand times 
the expected rate). 
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Figure 4: t + t~ qq channel: distributions of four analysis variables at different levels 
of the selection, as described in the text. Data at y/s = 191.6-201.7 GeV (dots) are 
compared with SM background expectations (left-hand side histograms) and with the 
expected distribution for a 105 GeV/c 2 signal in the (1i^t + t~)(Z— >qq ) (right-hand side 
histogram, normalised to one hundred times the expected rate). 
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Figure 5: HVz/ channel: distributions of four analysis variables, as described in the 
text, at preselection level. Data at y/s = 191.6-201.7 GeV (dots) are compared with SM 
background expectations (left-hand side histograms) and with the expected distribution 
for a 105 GeV/c 2 signal (right-hand side histogram, normalised to 10 times the expected 
rate). 
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Figure 6: Hz/P channel: Top: distributions of the likelihood variable for the expected SM 
backgrounds (full histograms), 191.6-201.7 GeV data (dots) and the expected Higgs signal 
at 105 GeV/c 2 (dashed histogram, normalised to 100 times the expected rate). Bottom: 
curve of the expected SM background rate at yfs = 199.6 GeV as a function of the 
efficiency for a 105 GeV/c 2 Higgs signal when varying the cut on the likelihood variable. 
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Figure 7: Hqq channel: distributions of four analysis variables, as described in the text, 
at preselection level. Data at y/s = 191.6-201.7 GeV (dots) are compared with SM 
background expectations (left-hand side histograms) and with the expected distribution 
for a 105 GeV/c 2 signal (right-hand side histogram, normalised to the expected rate). 
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Figure 8: Hqq channel: Top: distributions of the neural network variable for the expected 
SM backgrounds (full histograms), 191.6-201.7 GeV data (dots) and the expected Higgs 
signal at 105 GeV/c 2 (dashed histogram, normalised to 20 times the expected rate). 
Bottom: curve of the expected SM background rate at y/s = 199.6 GeV as a function 
of the efficiency for a 105 GeV/c 2 Higgs signal when varying the cut on the neural 
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Figure 9: hA hadronic channel: distributions of four analysis variables, as described 
in the text, at preselection level. Data at y/s = 191.6-201.7 GeV (dots) are compared 
with SM background expectations (left-hand side histograms) and with the expected 
distribution for a 85 GeV/c 2 signal at tan/3 = 20 (right-hand side histogram, normalised 
to the expected rate). 
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Figure 10: hA hadronic channel: Top: distributions of the likelihood variable for the ex- 
pected SM backgrounds (full histograms), 191.6-201.7 GeV data (dots) and the expected 
85 GeV/c 2 Higgs signal at tan/3 = 20 (dashed histogram, normalised to 20 times the 
expected rate). Bottom: curve of the expected SM background rate at yfs = 199.6 GeV 
as a function of the efficiency for a 85 GeV/c 2 Higgs signal at tan/3 = 20 when varying 
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Figure 11: Distribution of the reconstructed Higgs boson mass(es) when combining all 
HZ or hA analyses at 191.6-201.7 GeV. Data (dots) are compared with SM background 
expectations (full histograms) and with the normalised signal spectrum added to the 
background contributions (dashed histogram). Mass hypotheses for the simulated signal 
spectra are mn = 105 GeV/c 2 in the HZ channel and tba = 85 GeV/c 2 , tan/? = 20 in 
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Figure 12: SM Higgs boson: confidence levels as a function of mn . Curves are the 
observed (solid) and expected median (dashed) confidences from background-only exper- 
iments while the bands correspond to the 68.3% and 95.0% confidence intervals from 
background-only experiments. Top: 1- CLb for the background hypothesis. Also shown 
here is the curve of the median confidence as expected for a signal of mass given in the 
abscissa (dotted line). The sensitivity for a 5a discovery, defined by the horizontal line 
at 5.7 10~ 7 , is for Higgs masses up to 98.2 GeV/c 2 . Bottom: CL S , the pseudo-confidence 
level for the signal hypothesis. The intersections of the curves with the horizontal line at 
5% define the expected and observed 95% CL lower limits on mn . 
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Figure 13: SM Higgs boson: test-statistic Q for each mn hypothesis in data (solid) and 
its expected median value in background-only experiments (dashed). Top: the bands cor- 
respond to the 68.3% and 95.0% confidence intervals from background-only experiments. 
Bottom: the bands represent the confidence intervals around the minima of the —2lnQ 
curves expected for a signal of mass given in the abscissa while the dotted lines show the 
median curves expected for three mass hypotheses. 
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Figure 14: 95% CL upper bound on £ 2 , where £ is the HVV (V=W ± or Z) coupling 
normalised to that in the SM, assuming SM branching fractions for the Higgs boson. 
The limit observed in data (solid line) is shown together with the expected median limit 
in background-only experiments (dashed line). The bands correspond to the 68.3% and 
95.0% confidence intervals from background-only experiments. Most of the jumps below 
40 GeV/c 2 correspond to the different analyses applied to different subsets of LEP1 
data to cover the various topologies expected from a Higgs boson. The jumps at 12, 50 
and 65 GeV/c 2 indicate where the LEP2 results at 191.6-201.7, 182.7 and 188.7 GeV/c 2 
respectively start to contribute. 
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Figure 15: hA analyses: test-statistic (top) and confidence levels in background hy- 
pothesis (middle) and in signal hypothesis (bottom) as functions of m^+rriA- Curves are 
the observed (solid) and median expected (dashed) confidences from background-only ex- 
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Figure 16: MSSM Higgs bosons: regions excluded at 95% CL by the searches in the hZ 
and hA channels up to ^/s = 201.7 GeV, in the m™ ax scenario. The dark shaded areas 
are the regions not allowed by the MSSM model in this scenario. The dashed lines show 
the median expected limits. 
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Figure 17: MSSM Higgs bosons: regions excluded at 95% CL by the searches in the hZ 
and hA channels up to ^/s = 201.7 GeV, in the no mixing scenario. There is a region 
at mh around of 85 GeV/c 2 and small tan/5 that is not excluded, but is too small to be 
visible in the top left-hand plot. The dark shaded areas are the regions not allowed by 
the MSSM model in this scenario. The dashed lines show the median expected limits. 
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Figure 18: MSSM Higgs bosons: regions excluded at 95% CL, by the searches in the 
hZ and hA channels up to \fs = 201.7 GeV, in the large (i scenario. The dark shaded 
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Figure 19: MSSM Higgs bosons: regions excluded at 95% CL by the searches in the hZ 
and hA channels up to y/s = 201.7 GeV, in an extended scan of the MSSM parameter 
space. The dark shaded area is the region at high m h not allowed by the MSSM model 
in this scan. 



